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ABSTRACT 



Context. Recently, the MAGIC collaboration reported evidence for a delay in the arrival times of photons of different energies during 

a y-ray flare from the blazar Markarian 501 on 2005 July 9. 

Aims. We aim at describing the observed delayed high-energy emission. 

Methods. We apply a homogeneous synchrotron self-Compton (SSC) model under the assumption that the blob containing relativistic 
electrons was observed in its acceleration phase. 

Results. This modified SSC model predicts the appearance of a y-ray flare first at lower energies and subsequently at higher energies. 
Conclusions. Based on the reported time delay of ~ 240 s between the flare observed at ~ 190 GeV and 2.7 TeV, we predict a delay 
on the order of 1 h if observed between 10 GeV and 100 GeV. Such delay timescales can be tested in the future by simultaneous flare 
observations with the Gamma Ray Large Area Space Telescope (GLAST) and Cherenkov telescopes. 
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1. Introduction 

Recent observations of energy-dependent, few-min ute timescale 
TeV y- ray flares in the blazar Markarian (Mrk) 501 (lAlbert et al.l 
l2007al) are diflicult to explain in terms of the "classical" 
homogeneous synchrotron self-Compton ( SSC) model (e.g., 
iMarscher & Geadl985tlMaraschi et all 1992) in which the emis- 
sion regio n is moving along th e jet with a constant Lorentz factor 
~ 10-20 (iGiroletti et al.l2004 . as inferred from the observations 
of superluminal motion in active galactic nuclei. Similar few- 
minute timesca le flares have also been observed recently from 
PKS 2155-304 dAharonian et al.ll2007h . However, in the case of 
PKS 2155-304, no energy-dependent flare evolution was found. 

Homogeneous blob models will require blobs moving with 
Lorentz factor s significantly larger than m entioned above, see, 
e.g.. Fig. 3c in Bednarek & Protheroe ( 1999), or the more recent 
papers by Begelman et al. (2008) and Aharonian et al. (2007). 
In the case of the 2005 July 9 flare of Mrk 501, an arrival time 
delay of 239 ±78 s was observed between y-rays in the > 1 .2 TeV 
energy range with r espect to the 0.15 - 0.25 TeV energy range 
dAIbertet al.ll2007al) . This delay was confirmed by a more de- 
tailed quantitati ve investigation wit h an unbinned, photon-by- 
photon analysis d Albert etal.ll2007bl) . 

Two solutions for these fine emission features have been con- 
sidered u p to now, viz. gradu al acceleration of electrons inside 
the blob (lAlbert et aT]|2007al) and the attri bution of the dela y to 
a source-extrinsic y-ray transport effect (Alber t et al.ll200 7?). In 
the gradual-acceleration scenario, particles inside the emission 
region moving with constant Doppler factor need some time to 
be accelerated to energies that allow them to produce y-rays 
with specific energies. Therefore, a y-ray flare should first ap- 
pear at lower energies and only later on at larger energies. The 
second scenario relates the delay between different energies to 
the energy-dependent speed of photons in vacuum. The effect 
that lower-energy photons aiTive before the higher-energy ones 



is pred icted in some models of quantum-gravity (see lAlbert et al.l 
I2007bi for details). 

The large Doppler factors inferred from TeV y-ray observa- 
tions have to be achieved somewhere in the very inner parts of 
the jet. Therefore, it is natural to expect that the emission region 
has to undergo an initial acceleration phase. Very -Large Baseline 
Interferometry (VLBI) observations ra ther show evidence for 
decelerating radio-emitting blobs (e.g., lEdwards & Pinerll2002l 
for Mrk 501). However, the very inner regions where we assume 
the initial accereration to take place, could not be investigated 
yet with VLBI. 

We show that a delayed emission of the higher-energy y-rays 
can be naturally explained by assuming that the blob acceler- 
ates significantly during its propagation along the jet, produc- 
ing more energetic y-rays when the Lorentz factor of the blob is 
larger. We argue that for the first time y-ray emission in BL Lac 
objects has been captured in the initial acceleration phase of the 
blob. 



2. Accelerating blob model 

We consider a blob of relativistic electrons that moves along the 
jet with a Lorentz factor yb at any specific moment. We assume a 
"classical" SSC model, in which electrons are accelerated inside 
the blob by shock acceleration, with an acceleration efficiency 
proportional to the magnetic field strength at the location of the 
blob. Due to the presence of magnetic fields, the electrons inside 
the blob lose energy by synchrotron emission. This soft radia- 
tion is upscattered to y-ray energies by the same population of 
electrons. In terms of this classical picture. X-ray and TeV y- 
ray flares are naturally explained. However, the observation of a 
time delay between y-rays of different energies belonging to the 
same flare cannot be easily understood in such a simple scenario. 
Therefore, we modify the homogeneous SSC model by assum- 
ing that the blob has been captured during its acceleration phase 
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in the inner part of the jet. Let us consider the blob accelera- 
tion between its Lorentz factors y™" and y™^". We assume that 
the magnetic field strength in the jet (and also in the blob) drops 
proportionally to the distance from the base of the jet, B cc X, 
which is a good approximation in the case of a conical jet with 
a fixed opening angle at some distance to the base of the jet. In 
our model, we discuss a quasi-spherical blob with a radius 



R^ocX, 



(1) 



i.e., proportional to its distance from the base of the jet. This is 
expected for a conical jet, whose borders limit the expansion of 
the blob. 

We propose that a delay observed in a flare between different 
energy ranges is due to the fact that the lower-energy part of the 
flare is produced when the blob is closer to the base of the jet 
(moving with a lower Lorentz factor), while the higher-energy 
part is produced at larger distances from the jet base, at which 
the blob has accelerated to a higher Lorentz factor 

We now derive a description of the acceleration of the blob 
along the jet. Observations show that the TeV y-ray luminosities 
found for Mrk 501 during the 2005 July 9 flare i n the above- 
mentio ned two energy ranges are comparable (lAlbert et alj 
l2007ah . The observed y-ray luminosity is related to the lumi- 
nosity in the blob frame by 



^obs ^b' 



(2) 



where D\, is the Doppler factor of the blob. In the simplest pos- 
sible case, in which the observer is located on the axis of the 
jet, the blob Doppler factor is related to its Lorentz factor by 
Db = [yb(l -ySbCosfl)]"' ^ 2yb, where cosO = 1 is the cosine 
of the observation angle, assumed - 0°, and /3b = v'b/c is the 
velocity of the blob. Moreover, the y-ray luminosity in the blob 
frame depends on the energy density of the synchrotron radiation 
and the Lorentz factor of the electrons. 



(3) 



where ye is the Lorentz factor of the electrons that produce y- 
rays (with characteristic energies for which a time delay has been 
reported), and A^e is the number of relativistic electrons inside 
the blob, which is assumed to be independent of time. We calcu- 
late the average energy density of synchrotron photons inside the 
blob by estimating the synchrotron energy losses of the complete 
electron population inside the blob. That is, the energy losses of 
all electrons inside the blob multiplied by the average time that 
the synchrotron photons spend inside the blob, ~ Rt/c, divided 
by the volume of the blob. 



Rb 



yn 



where 



Vb 



Psyn ^ BY^N, 



(4) 



(5) 



are the synchrotron energy losses of electrons, Rb is the radius 
of the blob, and Vb oc R^ is the volume of the quasi-spherical 
blob. Provided that the energies of the electrons inside the blob 
are limited by their synchrotron energy losses, we compare their 
synchrotron energy losses with their energy gains from the ac- 
celeration mechanism, which scales 2LS -Pace ^ 

B in the shock ac- 
celeration scenario. Then the maximum energy of the electrons 
is 



In this case the synchrotron energy losses of the most energetic 
electrons are only linearly proportional to the magnetic field 
(Eq. |5]). The blob is assumed to expand during its propagation 
along the jet proportionally to the distance traversed, i.e. its vol- 
ume increases like Vb oc X^. From these simple considerations 
we conclude that the y-ray luminosity in the blob frame depends 
on the distance from the base of the jet as 



(7) 



Since the observed y-ray luminosities (P^^) in the considered 
energy range are approximately constant, the Doppler factor of 
the blob has to increase during the propagation of the blob along 
the jet like Db oc R^, to keep the y-ray luminosity approximately 
constant during the flare (cf. Eq.|2]l. Since the blob moves with 
a large Lorentz factor, the traversed distance is given approxi- 
mately by X ^ cf oc /^b, i.e., its Lorentz factor has to be approxi- 
mately proportional to the time during which it travels inside the 
jet. Note that the time t is measured in the stationary jet frame. 
Thus, 



yb ^ 0.5Db = At, 



(8) 



where A is a constant. We conclude that the observed time de- 
lay between different energies can be related to the accelera- 
tion phase of the blob inside the jet. In the next section, this 
conclusion is further investigated to constrain the blob accel- 
erati on scenario for the case of the recent observations of Mrk 
501 (lAlbert et alJl2007ah . 



3. Energy-delayed y-ray flare 

In our model, electrons are accelerated to maximum energies 
that change during the propagation of the blob due to a change 
of the local magnetic field strength in the jet. Thus, the y-rays 
produced by these electrons also might have characteristic ener- 
gies that are related to the energies of their parent electrons. The 
energy of a produced y-ray photon can be estimated by 



(9) 



This estimate is reasonable in the case of a relatively flat power- 
law spectrum of synchrotron photons (as observed during the 
Mrk 501 flares), i.e., when the inverse Compton scattering 
mainly occurs in the transition region between the Thomson and 
the Klein-Nishina regime. The time delay between the appear- 
ance of the flare at two different energies, is™" and is™", is due 
to the production of the corresponding y-rays at different stages 
of the acceleration of the blob. While a photon with energy /i™" 
moves the distance c df, the blob traverses a distance yScdf, where 
dt is the time measured in the jet frame, not in the blob frame. 
Therefore, the distance (in the jet frame) between two photons 
Ef" and E™"" is d{AX) = c(l -/3)dt. The distance between these 
photons corresponds to their time delay measured in the observer 
frame d(Ar) = d(AX)/c. The total time delay in the case of a 
blob with changing velocity can be calculated as 



At: 



(10) 



Je oc B 



-1/2 



(6) 



where 1 - pit) = l/2y^, f is the propagation time of the blob 
inside the jet, and f,ni„ and f,^ax denote the moments at which the 
Lorentz factor of the blob is equal to y™" and y™", respectively. 
By applying the relation between the Lorentz factor of the blob 
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and the travel time (Eq.[8j, we arrive at a formula for the obser- 
vational time delay between y-rays of different energies 



Ar 



1 rrr 
^ 2A 4„ 



dyb 



1 

2A 



2Ay\ 



mm 



(11) 



provided that y™^" » y™". This simple relation for time de- 
lays in the flare observed at different photon energies holds true 
as long as these energies are directly linked to the energies of 
the parent electrons accelerated in the blob. It is obvious that 
the lower the energy ranges are chosen, the larger the observed 
time delay will be: For example, a time delay between y-rays of 
30 Ge V and 300 GeV should be larger by an order of magnitude 
than a delay observed between y-rays of 300 GeV and 3 TeV. 
This clear prediction can be tested in the near future by simulta- 
neous observations in the GeV energy range (with GLAST) and 
the TeV energy range (with MAGIC, H.E.S.S., and VERITAS). 

During its acceleration from y™'" to y™", the blob covers a 
distance 



■ r"'"/J(f)df- ^ p (l-0.5yb2)dyb 



'b 



max 
'b 



7( 



Tb -Tb 



(12) 



in the jet. By using Eq.[TTl we can express the distance Xacc trav- 
eled by the blob by the measured time delay At and the limiting 
values of the Lorentz factors of the blob, 

X,,, = cAT(2y->™"-l). (13) 



4. Constraints on the Doppler factor of the blob 

Since a homogeneous SSC model is considered, we can con- 
strain the Doppler factor of the blob in two ways based (a) on 
the known variability timescales of the synchrotron and y-ray 
emission; and (b) on the escape condition of y-rays from the 
synchrotron radiation inside the blob. The maximum energy of 
the observed synchrotron photons can be estimated by 



B 



B 



Ocr Ocr 



rrnax 



D7 



(14) 



where B„ - 4.4 x 10'^ G and y^ can be approximated by apply- 
ing Eq. |9] ye ~ £'™''/(meZ)b)- On the other hand, the observed 
emission variability timescale of the flare has to be at least equal 
to (or larger than) the synchrotron cooling time of the electrons 
to observe a synchrotron flare, i.e.. 



Db' 



(15) 



where = lim^y^lCB^yl, and C is a constant describing the 
efficiency of synchrotron energy losses. By estimating the mag- 
netic field strength from Eq.[T4]and reversing Eq.[T5] we obtain 
a lower limit on the Doppler factor of the blob. 



(16) 



By considering the most extreme pa rameters ever observed 

20 TeV (lKonopelkoetal.l[T99l 



lAharonianet"ani2001h, s, ^ 5 MeV (ICatanese et al.lll997l) . 
and Tvai ~ 200 s dAlbert et al.ll2007a l). we estimate the lower 
limit on the Doppler factor of the blob to be ~ 10, and the 
lower limit on the magnetic field strength inside the blob region 
to be B Si 0.3 G. Note that these values are consistent with esti- 
mates deri ved from more complete info rmation on a flare in Mrk 
501 (e.g., iBednarek & Prothero3ll999 ) and with Doppler fac- 
tors ~ 25 - 50, which have been employed (Albert et al. 2007^ 
for the modeling of the Mrk 501 multi-wavelength spectrum 
for the MAGIC 2005 observations. Note, however, that model- 
ing the TeV y-ray emission from PKS 2155-304, in which also 
very short timescale flares have been obse rved, seems to require 
Doppler factors on the order of ~ 100 ( Aharo nian et al.l l2007). 

We can also infer a lower limit on the Doppler factor of the 
blob by requiring the optical depth for y-rays in the synchrotron 
radiation of the blob to be below unity. For another lar ge flare ob- 
served from Mrk 501, Bednarek & Protherod (Il999l) estimated 
the optical depth for y-rays produced in the blob. Based on their 
Eq. (6), using the above mentioned parameters for the 2005 July 
9 flare, we estimate the Doppler factor of the blob to be > 25. 
This rough estimate indicates that the blob has to move, at least 
during a part of its path inside the jet, with Lorentz factors above 
ten. Below we try to constrain the relation between the Lorentz 
factors of the blob during its acceleration phase by using the 
available information from observations in the context of the 
SSC model considered. 

In the discussion above it was assumed that the maximum 
energies of the accelerated electrons are related to the magnetic 
field strength at the blob location (see Eq. |6]l, because the syn- 
chrotron process constitutes the dominant energy loss mecha- 
nism. But they are also related to the maximum energies of the 
y-ray photons produced at the respective specific location of the 
blob inside the jet (see Eq.|9]l. In the considered model, the mag- 
netic field in the blob drops like B cc X ' oc R^^, with the 

Doppler factor of the blob D\, oc R^, and y^ oc B Therefore, 
the energies of the observed y-ray photons emitted at a specific 
location of the blob should be related to the Doppler factor of 
the blob by (see Eq.^ 



Ey OC 



,3/2 



(17) 



for Mrk 501, i.e., E 



y.max 



The MAGIC collaboration reported the observation of a time 
delay between y-rays at different energies /i™" and /i™^". For the 
observed flare of Mrk 501 on 2005 July 9, these limiting values 
differ by approximately one order of magnitude. Therefore, the 
following relation should be fulfilled by the Doppler factors of 
the blob at the two characteristic blob locations in the jet, which 
correspond to the production of y-rays with the energies Zi™" 

and E™": D^V^^r" ~ y^Vy"'" ~ 10^^^ 

As an example, based on the above constraints, we consider 
a blob that is in its acceleration phase characterized by Lorentz 
factors y"'" = 10 up to y"^" - 50. For the reported time delay 
between the flare at different energies (Ar = 239 s), we can esti- 
mate (from Eq.[TT]i the value of the parameter A ^ 1.7xl0"^s"'. 
The distance traveled by the blob during the acceleration phase 
between the two considered locations in the jet is Zacc ~ 7 x 10'^ 
cm (cf. Eq. [T2] |. Note that this distance scale is ~ 20 times 
larger than the Schwarzschild radius of the IO^Mq black hole 
(iFalomo et al.ll2002t iBarfli et alj|2003h expected in the center of 
Mrk 501. Therefore, in our model, the acceleration region of the 
blob containing relativistic electrons might well be located rela- 
tively close to the central engine. At such distances even the ra- 
diation field from the accretion disk around th e black hole might 
prevent the escape of y-ray photons (see, e.g.. lBecker & KafatosI 
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Il995h . However, any conclusion on such effects is limited by the 
unknown radiation fields provided by the accretion disk and the 
broad fine region in the case of BL Lac objects. 

5. Conclusions 

We have modified the classical homogeneous SSC model by al- 
lowing the blob to accelerate during its propagation along the 
inner part of the jet. The maximum energies of the electrons re- 
sponsible for the y-ray production are determined by the syn- 
chrotron process whose efficiency depends on the local strength 
of the magnetic field in the blob. Therefore, the maximum en- 
ergies of y-rays are also related to the acceleration phase of the 
blob. In the resulting y-ray flare, photons with lower energies 
should be observed prior to those with higher energies. Such a 
modified SSC model can explain the recently observed arrival 
time difference of y-rays belonging to different energies in a 
flare of Mrk 501. In our model, for the first time the MAGIC 
telescope was able to observe the acceleration phase of a blob in 
the jet of Mrk 501 . Note that currently other mod els of the emis- 
sion features at VHE y-rays dAlbert et al.ll20073lbl) also provide 
consistent explanations. 

Based on the measured time delay between the flares at dif- 
ferent energies, we estimate the distance scale on which the ac- 
celeration of the blob occurred and predict that the time delay 
should increase inversely proportional to the energy of the y-ray 
photons. For example, the delay between the flare observed be- 
tween 20 Ge V and 200 GeV should be approximately an order of 
magnitude larger than the observed delay between 200 GeV and 
2 TeV. This clear prediction can be tested by simultaneous ob- 
servations of strong flares at GeV and TeV y-ray energies using, 
e.g., the LAT instrument on board of GLAST (at GeV energies) 
and Cherenkov telescopes (at TeV energies). For the specific 
case of the observed Mrk 501 flare, where a ~ 240 s time delay 
has been reported between ~ 190 GeV and 2.7 TeV (weighted 
average energies of the reported energy ranges with a power-law 
slope of a = -2.2), we predict a corresponding ~ 0.5 - 1 h de- 
lay between the peak position of this flare observed at 20 GeV 
and 200 GeV. Observations of such delay timescales are quanti- 
tatively consistent with the picture of an accelerating blob within 
the jet. 
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